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Abstract 



J^j ■ In the paper we consider non-relativistic analog of the Nambu-Jona-Lasinio 

{Sj . model and demonstrate, by this most simplified example, the method of dynam- 

ical mappings of Heisenberg fields on "physical" fields. We obtain the expressions 
for one particle eigen-excitation spectrum of full Hamiltonian (spectrum of "clad" 
fermions), wave function of two excitations' bound state and its mass, and compute 
the Bethe-Solpither amplitude. We derive as well the relation for the density of 
vacuum energy as a function of mapping parameters. The minimization of energy 
leads us to the fixation of these parameters. The method admits the generalization 

Qh! on the relativistic case. 

CD 



Introduction 

The present paper is devoted to the solution of the Heisenberg equation for four- 
fermion contact interaction (non-relativistic analog of the Nambu-Jona-Lasinio (NJL) 
model without form factor |1], Q]) in representation of physical particles and their states. 
The idea of "physical" states was first introduced by Heisenberg 0. However, the most 
rigorous and consistent realization of the physical states is given in the papers Q]. The 
Hamiltonian of the interacting particles system on these states by definition has the 
diagonal form: 

(a\H\b}=<a\ J d 3 qA + (q)A(q)E(q) | b > +W < a \ b > . (1) 

The important fact is that the states A + (k) | 0) are eigenstates of the Hamiltonian H, 
that is 

~H,A + (k)] |0) = E(k)A + (k)\0) (2) 

A(k)\0) = 0, 

where E(k) determines thus the spectrum of one particle excitations of "physical" vacuum, 
defined by the second line in (|[). 

Since all interaction characteristics are reduced to the calculating of matrix elements 
of Heisenberg field combinations over the "physical" vacuum, it is necessery to have a 
connection between the Heisenberg fields and the "physical" ones, whose excitations are 
described by the A + (k) operators. We will call such a connection, following the papers|Q, 
the dynamical mapping of Heisenberg fields on "physical" fields. The regular procedure 
of seeking the dynamical mappings is based on the fact that they must be consistent both 
with the Heisenberg equation for interacting fields and with relation (0). The meaning 
of the latters, in fact, is led to the statement that the state ^4 + (k) | 0) describes free 
quasiparticle with the momentum k and the energy E(k). It is worth noting that the 
mehtod of the dynamical mappings is naturally adjusted to describing the phase transition 
phenomena and the processes with the spontaneous symmetry breaking. This fact is 
connected with the known case, when the dynamical mappings correspond to the unitary 
non-equivalent transformations of the Hilbert space, transfering it to the space ortogonal 
to the initial one, accompanied by the reconstruction of the Hamiltonian spectrum and, 
in particular, by changing the vacuum (ground state) energy. We show below that the 
vacuum energy density depends on the parameters contained in the dynamical mapping. 
By varying the values of these parameters we " enumerate" all possible mutually ortogonal 
Hilbert spaces. One of these spaces is the "physical" space, and the criterion of its choosing 
is the condition for the vacuum state energy to be minimal. 

Further on we consider the bound state of two excitations (of two "clad" fermions). 
It is well known, that while investigating the bound states properties in quantum field 
theory the main quantity is the Bethe-Solpither (BS) amplitude |5|], given by the following 
matrix element: 

$ afec (x, y; k) = (0 | TA a (x)B b (y) | A c (k)} , (3) 



where A a (x) and Bb(y) are field operators in the Heisenberg representation, and | A c (k)) 
is the bound state vector, a,b,c- generalized indices (spin, flavour, colour, ect), depending 
on the choice of studied theory. The usual approach is to use the BS equation |J ether on 
®abc(x, V] k) or on vertex function: 

T ab ...{x 1 ,x 2 , . . .) = (0 | TA a (x 1 )B b (x 2 ) ... | 0) . (4) 

The type of the vertex is fixed by the interaction. It is worth mentioning, that reactions 
of composite particles can be described by the effective Lagrangians (see, for example, 
reviews |], ?]). Development of QCD initiated the interest to the models with four-fermion 
interaction || (non-local prototypes of the NJL model) in which gluons are effectively ac- 
counted by introducing the non-local quark interaction kernel [§]. 

The equations on the function (H) or (||) are usually solved by using some perturbative 
theory (loop expansion, "ladder" approximation, ect). This leaves the set of crucial points 
out of consideration. Firstly, the matter of vacuum state choice. Secondly, connection 
of initial field operators in Heisenberg representation with operator A + (k), forming the 
bound state. At last, for the solution of BS equation it is necessery to know total two- 
point Green functions, i.e. the solution of the Schwinger- Dyson (SD) equation. The 
method of dynamical mappings allows to overcome all these dificulties, because it gives 
the opportunity to do direct calculations of all interaction physical characteristics. 



1. Dynamical mapping of the Heisenberg fields 

The Hamiltonian of the four-fermion interaction can be written in the following form: 



H = d 6 x 



V£ (x)e(x/)^ a (x) + -X + (x)x(x) 



(5) 



where a is spin index running over 1,2, e(v) " energy spectrum of free fermions, defined 
by the condition 

e ( v )e ikx = e(k)e ikx (6) 

and 

X(x) = e a 0ip a (x)ip p (x), x + (x) = e a(S i/j^(x)ip+(x). 

We need to find the relation between the Heisenberg fields i/) a (x) and fields (f> a (x), whose 
excitations would be eigenstates of the Hamiltonian H. In other words, it is necessary to 
find out the dynamical mapping [|J. As a starting point we consider the following case. 
Let the fields ip a (x) satisfy the canonical equal-time anticommutation relations: 

{tPa(x),^l3(y)} tx=ty =8ap8(x-y). (7) 

Defining generator W as: 

W= fd 3 x (f(x) X (x) - 7{x)x + (x)) , (8) 



where f(x) is arbitrary complex function, we can write down the transformation of the 

fields 1p a {x) — > <fia{x) 

(j> a {x) = e w ij a {x)e~ w . (9) 

In this equation all the fields are taken at equal times. Using equation fl7|), we obtain: 

<j> a {x) = u(x)ip a (x) - v(x)e a 0ip^(x) 

if> a (x) = u{x)(f) a (x) + v(x)e a p4>^{x) (10) 



where 



f( x ) 
u(x) = cos | 2f(x) |, v(x) = sin \ 2f(x) 



From the explicit form of u(x) and v(x) there follows the normalization condition: 

u 2 (x)+ | v 2 (x) |= 1 (11) 

Second equation in ( [10f) gives rise to the dynamical mapping of the Heisenberg fields i/} a (x) 
on the physical fields (j) a {x) if, certainly, it is consistent with the Hamiltonian @. By the 
next step we postulate the set of properties for the field (f> a (x). That is, let 



[x 



(27T) 



■f d 3 kg(k)e lkx - iE{k)t A a (k). (12) 



be a representation for 4> a (x), where ^(k^nd A+(k)are annihilation and creation oper- 
ators for the particles with momentum k, satisfying the following canonical relations: 

{A a (k),A$(q)} = 6 a0 6(k- q). (13) 

-E'(k) is unknown yet excitation spectrum of A^ (k); g(k) is a distribution function over 
the momenta inside the packet, and g(— k) = g(k). Here we should note the following 
important circumstance. Representation ( |T2"D completely breaks the initial properties 
of transformations @ and flip]) - they cease to be canonical, i.e. they do not preserve 
the anticommutation relations. It is the presence of g(k) function that causes the field 
anticommutation relations of <j) a (x)not to be local. Indeed, from equation (|12|) with 
account of relation (|I3"D we have: 

{*«(*), #(V)} = ^rff d3k I $( k ) I' eWwH^WM] 

{M*)A + M x = y = S -^ (14) 

{0 a (x),A+(k)} = -^ T ^(k)e ikx . 
L J (2%) 2 

where 

' ' r d 3 k\g(k) | 2 (15) 



V* (2tt) 3 

4 



So, the input of g(k) allows to modify the product of fields to be defined at one point, thus 
breaking the locality. Transformation (||) plays auxiliary role for obtaining fllOl) leading 
to the fact that relations (|12|) and (fEf ) have no sense any longer. 

To find the meaning of normalization (|l|) we can use the dimensions of the fields 
(as it follows from (p~3|)) L4 2 (k)] = L 3 , ["0 2 (x)] = L~ 3 , therefore, g(k) is dimensionless 
quantity, and integral in the l.h.s of (|l|) has the dimension L~ 3 . Further on, since 
| g(k) | 2 characterizes the momentum fluctuations inside the excitation, volume V* has 
to be connected with a characteristic spatial size of the excitation. If one accepts that in 
(AA;) 3 region, where | g(k) | 2 essentially differs from zero, it is equal to 1 (| g(k) | 2 ^ 1), 
then 

pip/AUMI^fA^^ (16) 

where L* = j^- is linear size of the packet. Accounting this fact we take for | g(k) | 2 the 
normalization (|15|) with the additional condition: 



fd 3 k\g(k) 
jd 3 k\g{k)\* 



<l 000 \ 2 > 9 = V M l.)Z, = 1- (17) 



Thus, in momentum space typical for the wave packet <f)^(x), | g(k) | 2 ~ 1. 

So, we would search a solution in the form of relation (|i0|) with conditions (|12D , (p!4|), 
(|T5| ) and(p!7|). Let us integrate the left and right parts of the equation (|T0| ) with the 
plain wave and obtain expressions of Heisenberg operators a a (k) via creation v4+(k)and 
annihilation ^4 Q (k)operators of quasi-particle: 

a a {k) = fd 3 k'g{k')u(k-k')A a (k') + fd 3 k'g(k')v(k- k')e a pA+(k'). (18) 

Hence we can see that transformations (0), in contrast to the Bogolubov transformations, 
are nonlocal over momentum. Therefore, vacuum state | > defined as ^^(k)! >= 0, 
should essentially differ from the ground state defined by the Bogolubov transformations. 
The question is, whether the transformations ( |T0| ) are "compatible" (consistent) with the 
coupling dynamics of the Hamiltonian (|^)? It will be shown later on that at a particular 
choice of u(x) and v (x) such a compatibility takes place. 



2. One-particle excitation spectrum 

Let us choose the parameters u(x) and v(x) as follows 

u(x) = u = const, v(x) = voe ta ( x \vo = const 

ul + vl = 1 (19) 

Then the dynamical mapping will be of the form: 

ip a (x) = u <p a (x)+v e ia( - x) e af3 (j)p(x) 

V£(x) = u <Pi(x) + v e- ta ^e af iMx) (20) 



Our aim is to make the state A+ (k) | 0) to be an eigenstate of the Hamiltonian ([|), with 
the vacuum being defined with respect to the operator A a (k) 



A a {k) | 0) = 
Acting by the Hamiltonian (|5|) to the one-particle excitation we have: 
i?A+(k) | 0) = A + a {k)H | 0) + k^(k)l | 0) . 



(21) 



f22) 



If the vacuum is the Hamiltonian's eigenstate, then for the one-particle excitation to be 
stationary it is necessary to satisfy the condition: 



H,At(k) | 0} = E(k)At(k) | 0) 



(23) 



where E(k) is the energy of this excitation above the vacuum | 0). This relation can be 
thought either as a definition of the excitation spectrum or as the Heisenberg's stationary 
equation. 

Leaving for the moment aside the problem of stationarity of the vacuum, let us consider 
only the equation (|23|). The commutator in the l.h.s. of (^) splits to the commutators 
with the kinetic term and the coupling term of the Hamiltonian independently, leading 
to the relations: 



X + (x) X (x), A+{k)\ = X + {x) [x(x),A+(k)\ - [A+(k),x + (x)\ X {x), 
2u 



x(x),A+(k) 
A + a (k), x + (x 



(2vr 



g{k)e l ^e^p{x) 



(24) 



2v(x) 

(2vr)§ 



g(k)e^+(x). 



To derive fl2~4] ) we have used the dynamical mapping (PU[). Putting together all three 
relations of (f24T) we obtain: 



X + (x) X (x),A+(k) 



+ 



2ttp 
(2tt)I 
2v(x) 
(2tt)I 



g(k)e lkx X + (x)e 0a M x ) + 



g(k)e^:(x)x(x). 



So, the problem is now reduced to the calculation of two terms in the r.h.s. of 
the first one we get the following: 

X + (x)ip J (x) | 0) = v(x)e 11 ix + {.x)(l)^(x) | 0) = 

= v{x)e^ l e a ^{x)ip^{x)(j)^ l {x) | 0) = 

= v(x)e ri '6 a (u (j)p(x) + w(x)e /3/3 /0 /3 /(a;)J (u (f)+(x) + 

+v{x)e aa t(f) a r(x)) <j)y{x) I 0) = 

= v(x)e^>e a/3 (u v(x)e aa '(f)^(x) {<^>(x),0+,(x)} + 

+M O TJ(x)e / 3 /3 /0 /3 '(x)0+(x) 0y(x)) | 0) = 
UqvI 

v* 



(25) 
For 



77'£ Q /3 U a y4>^(x) - e^/(/»+(x) + ep a (f)^{x)) | 0) = 



(26) 



since the last expression in brackets is fully antisymmetric by indices ap'j. Consider now 
the second term of (|25|). 

i + xW | 0) = e^ip+ix^ix^ix) | 0) = 

= -v(x)^(x)tpp(x)^(x) I 0) = -v(x) (uo(f>+(x) + v(x)e aa xj) a i(x)^ x 

x (u Q (t}fi{x)(j)^{x) +v(»)e / j j9 /0j(x)0j(x)) | 0) = 

= -v(x) \2^<pl(x)\ - (y^e aa >e 0i<l>^(x)(f> a >(x)(l>j{x)) | 0) = 
2v(x) 



V* 



-#(a0|o). 



Here we have used the constraint it 2 , + i> 2 = 1 With account of the results (|26] 
the expression Q25D we will have: 



-X + (x) X (x),Ai(k) 



0) 



. A 2 s( k ) „ lkx 



y*(2vr 



e* kx 0+(x)|O>. 



(27) 
and 

(28) 



Integrating this relation over the whole configuration space the contribution into the 
energy spectrum of the Hamiltonian reads 



A 



d 3 x X + (x)x(x),A+(k) 



0) = -A^ | <7(k) | 2 A+(k) | 0) 



(29) 



Now we derive the kinetic term contribution: 



V^Mv)Vv(x), A+(k)] | o) = -^(x)A+(k) e (v)Vv(*) I o) + 

V>+,(x)e(v){VvOr),A+(k)} | 0) -A+(k)^+(:z)e(v)^'(*) I 0) = 

= -{^+(x),A+(k)}e( V >Mx) I 0) +^,(x)6( V ){^(x),A+(k)} | 0) = 

= -u(x) {<^(x), A+(k)} e( V Mx)0+(x) | 0) + 

+ M 2 0+(x) e ( V ) {^(x),X+(k)} | 0) = (30) 

= -v(x)^ T e^e( V )v(x)<f ) i(x) | 0) + ^-^W kx e(k)«^(x) | 0) 

(27TJ 2 (27TJ 2 

The last term in (|30|) when being integrated over the space variables is proportional 
to the operator ^4+(k),and gives therefore normal c- value contribution into the energy 
spectrum. Rather complicated situation arises after integrating the first term. Evidently, 
the diagonalization of this term cannot be done with an arbitrary function v(x). So, 
we need to find the conditions under which the diagonalization is possible. Taking into 
account that v(x) = Voe ia ( x > we can write: 



v(x)e(sj)v{x)(j)+ (x) = Uo e (v)0a(:r) + v(x) [e(\j),v(x)} <j>+(x) 



(31) 



Here square brackets mean commutator we need to calculate. 



Let f(x) be an arbitrary function , then one may write 



— [e(s/),v(x)]f(x) 



<W),e 



ice(x) 



/(^) = -2^[(v) 3 ,e fa(8) ]/(x) 



2m 



2 y e ia{x) v /(&) + /(s)(v) 2 e ia(x) 



Further we have 
therefore 



[V) 2 e ia{x) = e ia{x) (*(v) 2 a(x) - (V«W) 2 ) • 



Mv),v(x)] 



V{X) 

2m 



2(ya(i))k - i y 2 a(x) + (va(i)) ; 



(32) 



where k = — z'y. 

From the considerations above it can be seen that a diagonalization is possible only 
at ya(x) = const. In this case the coefficient at 0^(x)will be homogeneous and the space 
integration leads to the term proportional to A+(k). Keeping in mind all this we take 
v(x) in the form 

a(x) = k • x (33) 

where k is some constant vector, and its physical meaning we will explain below. For 
the commutator we finally get 



[e(v),v(aO] 



v{x) 
2m 



2k k + k 2 ) 



(34) 



Now we are able to write out the kinetic term contribution into the energy spectrum 

^W<V)i'W,^(k)] I 0) = ^e(k)-^e* kx 0+(x) | 0) - 

' v2 ° e{k) wr/^ ( e(v) + i (2ko ■ k + k ° } ) ^ {x) ' 0) (35) 

2 

Integrating this relation along the whole space and accounting that 



e(v) + 7T-(2k • k + kg) ) e"* kx = e(k - k )e- kx 
2m ) 



(36) 



we have 



d 3 ^+(x)e(y)Vv(z),A+(k) 



0) 



(37) 



= |^(k)| 2 (^ e (k)-^ e (k-k ))^(k)|0) 
Combining fl3"7| ) and (2~D) we obtain the expression for the energy spectrum of excitation 

E(k) =\ g(k) | 2 (u%e(k) - v 2 e(k - k ) - X^j (38) 



Thus, transformation (|20| ) under condition ( j33|) realizes, indeed, the dynamical mapping 
in a sense that the fields 4> a {x), entering the mapping, describe the Hamiltonian ([D 
eigenstates. However, all this is true due to the fact that the state without excitations 
| 0) is the Hamiltonian eigenstate as well. Let us dwell on this question. 



3. Vacuum stability and conterterm. 

Let us compute the action of the Hamiltonian on the state | 0) without excitations 

A 



H | 0) = f d 3 x 
For the kinetic term we have 



V£ {x)e{\/)tp a {x) + '-x + (x)x(x) 



V£(zMvW«(aO I o) 



0) 



(39) 



= e Q/3 (u (j)t(x) + v(x)e a i3>(j)p(xy) e(\/)v(x)(f)^(x) | 0) = 
=v(x)(t) a {x)e(\j)v{x)(t)l{x) | 0) + u e a p(j)+(x)e(x7)v(x)(j)^(x) | 0) . 

The corresponding expression for the interaction term brings the following result 

X + (x)x(x) | 0) = -v(x)x + (x)ip a (x)(j)+(x) | 0) = 
= -v(x)x + (x) (u (j) a (x) + v(x)e a/3 (p^(x)^ 0+(x) 
2u v(x) 



(40) 



V* 



-X + (x) | 0) -v 2 (x)e al3 x + (x)(j)i(x)(j)+(x) 



0) = 
0) . 



(41) 



For the first term in this sum we have 



0} 
0) 



x +(x) | 0) = e a 0i/)$(x)i/>+(x 

al 3l{jp(x) (u (f>t(x) + V(x)e a/ 3'(j)f3'(x 

u e a f3 (uo(f)p + v(x)ep/3>(f)0>(xy) cj>t(x) | 0) 

2v(x)uq 



(42) 



= -ule a/3 (f)+(x)(j)^(x) | 0) - 
For the second one we obtain 



V* 



0} 



£af3X + (x)<Pp(x)<p+(x) | 0) = e a pe al p4^,(x)ip+,(x)4>^(x)4>l(x) | 0) 
= e a pe a '/)'ipp(x) (u (j)a'(x) + v(x)e a i a ii(j) a „(x)j 4>^{x)4>+{x) | 0) = 



(43) 



v(x)e a ^,(f)0 l (x)(j)^(x)(j)+(x) | 0) 



L fl 2 (*)i n v 2 



0) - — u v{x)e c 



x 



(x) I 0) . 



Putting together (|4]), (0) and (f43"D we get 



X + (x)x(x) | 0) 



4 < ' - 2 



V 



, , 0) + — UQv(x)e a p^{x)^{x) | 0) 



(44) 



Thus, for the action of the Hamiltonian density on the vacuum we would have 

U+(x)e(y)ip a (x) + -X + (x)x(x)j | 0) = 

= v{x)t/> a {x) Uv) + 2^r) «(*)#(*) I 0) + (45) 

+u o e a0 (f)+(x) f e(v) + ^77) "W^W I °) 

Integrating this relation over the whole space we finally derive 

H\0) = consto I 0) + AH (2) | 0) (46) 

where 

const = J d 3 xv(x) I <f> a (x), (e(y) + 7^77 ) «O*O0£O»O > (47) 

and AH (2) is two-particle contribution equal to 

AH{2) = J d 3 xu e a ^t(x) fe(y) + ^M v(z)tf(z) (48) 

The braces in (|47|) stand for the anticommutator. Computing now consto by passing 
in ([17]) to the momentum representation and taking into account that e(k) is quadratic 
function of the momentum we can see 

2 r « , J . , < k 2 > A \ 



consto = — / d 3 xv (x) e(y) + 



2m 2V* 



= ^ / cftre^ (e(y) + ^ + ^) e*- = (49) 

n 2 V ( ,. , <k 2 > A \ 

Here V - is space volume, and 

< k >= j^U(k)l' (50) 

AH {2) being expressed in terms of the creation operators v4+(k)of quasiparticles has the 
form 



AH (2) = uov o e a0 J rf 3 x0+(x) fe(y) + ^j 



e Jk0X ^(x) 



/ to(k)e Q/3 A+(^ - k)A+(^ + k) (51) 



10 



where 



D(k) = u v g(k - ^)g(k + ^) (e(k h < 



2 w\ 2 / i v 2 ' 2K* 



Thus, as it follows from (^6[), the non-excited state | 0) is not an eigenstate of the Hamil- 
tonian (||) and the AH (2) term is the source of the nonstationarity of | 0). This term 
describes correlated fermionic couple of excitations, moving with the momentum ko. The 
physical meaning of the relation (f46|) is that it points to the existence of energy exchange 
between the couple and the system of fermions. In order to account this exchange one 
must input into initial Hamiltonian the term describing the couple. So, in the relation 
PB| ) we transfer AH(2) to the left part, thus, redefining the Hamiltonian and taking as a 
physical Hamiltonian the quantity Hp equal to 

H p =: H - AH{2) : (52) 

where the normal ordering is referred to the Hiesenberg fields. If 

AH{2) =: AH{2) : +const, 

then 

H p = H - AH{2) + const. H p | 0) = H | 0) - AH{2) | 0) + const | 0) , 

H p | 0) = W | 0) , Wo = const + const. (53) 

Consequently, after substitution of the conterterm AH (2) into the Hamiltonian (|5]) the 
non-excited state | 0) becomes stationary, with the energy equal to Wo- Will this redefi- 
nition of the Hamiltonian change the one-particle excitation spectrum -E'(k)? To answer 
the question let us write AH (2) in the following form: 



AH{2) = u o e a0 J d z x(j>+{x) fe(y) + 



A 




2V* 

u e af3 I d?xv(x)cj) a (x) e(y) + t^TT ) ^0*0 ( 54 ) 



The second term, we have added to the sum (|>|) , does not contribute to the relation ( f46D , 
however, providing the hermiticity of AH (2). The contribution of AH (2) to the energy 
spectrum is defined by the commutator (^3|) . It reads: 



AH(2),A+(k)\ |0>~A*|0>=0 (55) 

Thus, the addition of the conterterm AH(2) to the initial Hamiltonian does not affect 
the spectrum. 

It is worth noting that the Hamiltonian of Hp type has been studied in the papers (0), 
where some exact results for a set of physical quantities have been obtained. Our case 
has two essential differences. Firstly, AH(2) is added to the initial Hamiltonian not with 
an arbitrary distribution function D(k) but in the form "provoked" by the Hamiltonian 
H itself. Secondly, AH(2) describes correlated pairs not of the Heisenberg fields, but of 

11 



the excitations A+(k). We will show below that just these excitations form a bound state 
which is the ground state of the interacting fermions. 

At first sight AH(2) seems to break global gauge invariance. However, this is no true. 
The fact is that coefficient at D(k), connected with the rotation parameters uq and v(x), 
is expressed in terms of abnormal vacuum averages of the Heisenberg fields. So, under the 
gauge transformations of the fields these parameters are transformed as well by the way 
by which AH(2) remains invariant. In other words, condensate of correlated fermionic 
pairs restores the guage symmetry of the Hamiltonian. The constant entering (^) and 
defining the energy of vacuum state | 0) we will calculate later on, but now we consider 
two particle excitations of fermionic system with the Hamiltonian H p . 



4. Wave function and Bound State Spectrum of Fermionic Excitations. 

Let us find how the Hamiltonian H p acts on a state composed of two excitations with 



momenta k and q 

H p A+(k)A+(q) 



0) = (W + E(k) + E(ci))Ai(k)A+(c L )\0} + 
{[H pi A+(k)],A+}\0) 



+ 



(56) 



If the excitations do not interact with each other, the last term in the sum ( p6|) should be 
equal to zero. Then this two-particle excitation will be an eigenstate of the Hamiltonian 
Hp with the energy equal to the sum of energies of both excitations. Now we will demon- 
strate that this term does not vanish, reducing instead to the none-diagonal two-particle 
form. In this case, however, it will turn out that there exists a correlated combination 
of excitations which is stationary and defines the bound states of these excitations. The 
problem now comes to the computation of the last term in ([56]). Earlier above we have 
calculated the action of the commutator in (^) on the | 0) state defining, thus, the en- 
ergy spectrum E(k). It is necessary here to take into consideration those terms of the 
commutator as well, that vanished previously on the state | 0). All contributions to the 
commutator rise now from the kinetic term, interaction and conterterm AH (2). We begin 
treatment with the kinetic term: 



^+(s)e(vWaO,A+(k) 

V 

(27T 



'•< r)^W kx e a/3 e(v)<M*) + Uo ^e^e(k)^(x) 



(57) 



From which follows: 



{[^ 7 + (a;)e(v)^ 7 (x),A+(k)],A+(q)} = 
= e a ^(x)^^( e (q)+e(k))e^) 



(5f 



Let us note that relations (p?\) and fl58|) are absolute, i.e. they are fulfilled on any state. 
The analogous contribution coming from conterterm reads 



-Atf(2),A+(k) 



uov(x)e 11 ' 



4>t(x) e(y) + 



A 
2V* 



y'(*MJ(k) 
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u Q v{x)-^w^e^ a 



+ Uqv(x) 



(27T) 

<7(k) J 



kx. 



(2tt; 



C/3c 



e(k) + 
'e(v) + 



A 

2K*_ 

A 
2F* 



and 



{[-AJf(2),X+(k)],Aj(q)} = 

#Mq)i(k+q)x 



x + 



,r 



e^Mo^l^J 



(2tt) 



e(q)+e(k) + — 



(59) 



(60) 



Thus, as it follows from the relations ([58) and ( |60"D the contributions into anticommutator 
from the kinetic term and conterterm are just "c"-numbers. These contributions, as we 
will see later, compensate each other so that the only term proportional to the coupling 
constant A is left. Now we will show that this term, in its turn, is canceled by the 
contribution coming from the interaction term j (x + ( x )x( x )) °f Hp. As a result the 
sought anticommutator acting to the vacuum | 0) will be proportional to the two-fermionic 
state. Now the role of the conterterm AH (2) is completely determined, it annihilates two- 
particle component in the | 0) state and vacuum component in the A^,(k)At(q) | 0) state. 
Let us find out the contribution of the four-fermionic interaction. For the Heisenberg 
fields entering the relation (E3) one can write: 



2u 2 
ip+(x)x(x) = ule^,(f>+{x)(py{x)(f) 7 >{x) - ^-v(x)(f)+{x) + 

2 
+ 2u 2 v(x)(j)l(x) (j>^{x)(j)p{x) - —u o v o e a0 (j)0{x) + 

+ 2u o v%e a p<l>0(x)<f>+(x)<p 1 (x) + vlv(x)e a pe^<pp(x)(j)+ (x)<f)t(x) 



(61) 



It is easy now to calculate the action of the corresponding anticommutator on to the 
vacuum. Omitting the simple calculations, amounting to the commutation of A+(k)with 



fields (|6T|) we have: 



{ij+(x)x(x),A+( q )}\0)=v(x) 



glgj 
(2tt)I 



= ?qx 



e a /3e 77 >(f)+{x)(f)t(x) | 0) 



(62) 



Here we have used the fact that 



0a 0j I °) = 2 e «0 fc n'#y( a O#7'( a O I °) 



(63) 



Further, 



{x + (x)^ a (x),A+(q)} = x + (x){A^(ci),i; a (x)}+[A+(ci),x + (x)]ij a (x) 



{^3(q),^0)} = u 



(2tt)I 



Q '«qx; 



Jaj3 



(64) 
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Unification of these relations gives 

"A 



4 X + (x)x(x),A+(k) 



.4(q) |o) 



= ^ov(x)e a0 i^e^* | 0) + (65) 

for the interaction term. As it can be seen from this result the first addendum in (|65|) 
really compensates the last addendum in (|H]) so that the anticommutator acting to the 
vacuum gives the only two-fermionic state. Integration over the whole space yields: 

tf p A+( qi )A+(q 2 ) | 0) = (W + E(q 1 )+E(q 2 ))A+(cii)A+(< h )\0) + 

+ A ^ ( ^ ( 3 q2) J fe^+^ x 0+(x)0+(:r) | 0) (66) 

Here we would like to note that the excitations of two fermions with a total spin 1 
(symmetrical over a and 0) do not interact with each other, because in this case the last 



term in the sum (pq) is equal to zero identically and the energy of the state presents just 
the sum of free energies of either excitation. Let us pass to the variables P = qi + q2, 
k — \ (Qi ~~ 02)- Then using the expression for the fields <f> a (x) via operators A a (k), we 
can rewrite the equation (B6l) in the following form: 



P P 

|-k) + S(| + k)JA+(|-k)Aj(| 
A ,p , , ,P , , r * ,P 



H p AU--k)A+(- + k)\0) 
Wo + E(^ - k) + e£ + k) ) A + a {^ - k)A+(^ + k) I 0) + (67) 



r);r) 3^Y " %( J + k ) / d3 19(^ - CL)g(- + q)A+(- - q)A+(- + q) | 0) , 

Define the wave packet A + (P) as: 

i+(P) = J G(k, P)e aP A+(^ ~ k)Aj(? + k)d 3 k (68) 

and choose the wave function G(k, P) in order to A + (P) | 0) be an eigenstate of the 
Hamiltonian H: 

H P A + (P) I 0) = (Wo + MP)) A + (P) I 0) (69) 

After the simple calculations we obtain from (E7|) and (|69"D: 



C(k P)->v g(f-k)(7(f + k) 

MP) - E(f - k) - £(f + k) 



where 

P . , ,P 

To " 



Jd 3 kG(k,P)g(--k)g(- + k). 



(2tt 
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This set of relations determines the wave function G(k, P) and the energy of the two 
excitations bound state /x(P), with P being momentum of the bound state as a whole, so 
that /u(P = 0) is the bound state mass. In order to obtain the equation on /x(P) let us 
substitute the manifest form of the function G(k, P) (|70|) into the expression for 70. We 
see that 70 is cancelled from the left and right parts and we have the following equation: 

A te | 9 (f-k)P| g (f + k)P 



(2xfJ ^(P) - B(E - k) - B(f + k) 

This equation is the necessary and sufficient condition for existence of two particle bound 
state, and it determines the bound state energy /i(P). It is worth noting that this kind 



of condition always appears when searching the solution of the BS equation ||10|| , what 
follows from its homogeneity, and equation (|7T|) corresponds to the zero determinant of 
the homogeneous system. We will analyse this equation below, but now let us turn to the 
calculation of the BS amplitude $ a ^(x, y; P). It contains all the information about bound 
state. As has been pointed out in Introduction, to find it one has to solve the BS equation, 
which kernel in its turn is determined by the solution of the SD equation. However, if 
the dynamical mapping is known, it is possible to calculate $ Qi g(a;, y; P) directly, since, 
by definition, it is a matrix element between two physical states for the product of the 
Heisenberg fields: 



<f> a/3 {x, y; P) = (0 I T^ a {x)^ p {y) \ A(P)) , 

I A(P)) = i+(P) I 0) . (72) 

We can construct four BS amplitudes from the fields ip a (x). Introducing the variables 
R = £±£, r = x — y after straightforward calculations we can write for these amplitudes: 



I TiP+(x)Mv) I ^( p )) = u voe- ik ^e aa >S^G a/p/ (x,y;P); 

0\Ti> a (x)^(y)\A(P)) = uoVoe-^rewSacfGatfyixMP); 

(0 I Tip a (x)ip^(y) I A(P)) = ul5 aa ,5{3p,G a/ p,(x,y;P); (73) 

0\Tif>+(x)tf(y)\A(P)) = v 2 e- 2lkoR e aa ,e m G a ,p,(x,y-P), 

where the vertex part G aj g(x, y; P) is defined by the matrix element: 

G af3 (x, y; P) = (0 | <t> a (x)<f> (y) | A(P)) = (74) 

= -— ^- e iRP e Q/3 J d 3 kG(k, P)g(- + k)g(- - k) e ikr e-^ ( ? +k)4 ^ iS( ?- k) ^. 

Consider this relation at equal times t x = t y = t. We have: 

G aP (x,y;P) \ tx=ty=t = -^e iRp -^ p )*7og(r,t;P), (75) 

(27TJ2 

a(v fP)~ l f d 3 k ! 9 ^ + k) |2 ' g{ ^ - k) |2 ^ Am-E^-M)-E^ +k )) 
g(r ' tlP) " (2^)1 J d %(P)-i?(f-k)-i?(f + k) e 
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The plain wave appeared in (fTSP describes a movement of the bound state as a whole 
with the momentum P and the energy /i(P), and form factor g(r,£;P) is related with 
the bound state internal structure. At t — 0, P = q(r, 0; 0) = ^(r) coinsides with 
Schredinger wave function in "x" representation. Let us consider it in details: 

(2-k)2 J fi - E(k) - E(-k) 

is thought to be the bound state wave function in the configuration space. The bound 
state has the structure of the pair strictly correlated over the momenta. Parameter 70 is 
determined by the normalization of distribution G(k, P). Due to the fact that A < the 
following relation 

/j - E(k) - E(-k) < (77) 

should be satisfied for all k, since only in this case the signs of the left and right sides of 
the equation (|H]) will be conformed and the function under the integral will not have the 
peculiar points in the whole integration region. 

Let us calculate now the wave function $(x), considering the case when e(k) = ^— 
and representing G{k) in the following form: 



gik) | 2 M 



(7* 



Eik) + E(-k) - /i k 2 + 5 2 
Then the nonequality ([H]) corresponds to the conditions: 

M > 0, 5 > (79) 

By the use of the manifest form of the energy spectrum (|38|) we obtain: 

m 
M = , >0 

1 - 2vl 



., o . 2mA 2 my* 



<K + -Y7T< + —^- > 0. (80) 



2v 2 V ° ° V* ° ' I g(k) 



The approximation can be made by putting | g(k) | 2 = 1 and the integral ( [76|) will remain 
still convergent. As far as | ^(k) | 2 plays the role of a cutting-off factor at the large 
momenta, this approximation will lead to distortion of $(x) only in the region x ~ 0. 
Having in mind all mentioned above we can see, that from equation ( p0|) follows: 

1 - 2v\ > (81) 

Evidently, E(0) defines the lower boundary of the energy spectrum E(k). The relations 
(|8l| ) represent the necessary conditions for existing of the bound state of the fermionic 
excitations. There are two possible different physical situations according to the sign of 
£'(0). In the case when E(0) < the state | 0) with the energy Wo is located in the 
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excitation spectrum, and the bound state is the lowest, being at least 222(0) distant from 
the lower boundary of spectrum. Because 22(k) at k large enough becomes positive, there 
exists such a momentum po for which 22(po) = 0. Thus the state 

|0,l>=a|0) + c a A+(p )|0), (82) 

where a and c a - any complex numbers, is a stationary state with the energy Wq. If 
22(0) > 0, then for the value of /j, three variants are possible. Firstly, E(0) < ji < 222(0). 
Then the bound state lies within the excitation spectrum and the | 0) state with the 
energy Wq will be the ground state with the energy gap between it and the excitation 
spectrum. Secondly, when < \i < 22(0) two levels system arises with the ground state 
equal to | 0). And lastly, if \i < 0, we will have again two levels system, though the 
bound state will be the ground state of two excitations with the energy equal to \x. In 
the particular case when 22(0) = the momentum ko is determined by the four-fermionic 
coupling constant A. From the expression flST]) we obtain: 



, 9 2mA 

k^ = -yr (83) 

For the wave function $(x) we have, up to the normalization: 

$(x) = —L, f d 3 k — , [ 9 ^ [ e ikx ~ M [ d 3 k-^—— (84) 

(2tt)I ■/ 22(k) + 22(-k) - /i J k 2 + 5 2 K ' 



Integrating it over we get: 

$(x) = M e~ r5 (85) 

r 

As a consequence of the approximation we see the singularity of fljj5|) at r = 0. In order 
to account correction to (J35l) let us introduce the screening parameter r : 

2?r 2 
$(x) =► $(x) = M e~ rS (86) 



From here we obtain 



r = M 



r + r 
2n 2 



$(0) 



Further, as it follows from (|84j) 



* (0) = 7^1 / ft 



(2n)iJ 22(k) + 22(-k) - /i 

dk n\ 7 i = ~^" ( 87 ) 



(2tt)3 J 22(k) + 22(-k) - fj, A 

where we have used the equation (|7Jf) . Finally we derive 

AM 



ro = -17 
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It is seen from the solution ( p6[ ) that quantity ■h determines the radius of the bound state. 
Note, at last, that the wave function $(x) satisfies the equation 



(EM + E(- V ) - /i) $(x) = i(27r) 3 {^(x), 0+(O)} , 



(89) 



which follows from the manifest form of $(x). 

As a next step let us examine the solutions of the equation ([TI]) . Multiplication of the 
numerator and the denominator under the integral on the same quantity E(k) + E(— k) — /i 
gives for fi: 

H = 2(E(k)) a + W L^ (90) 



where 



2(E(k)) 



XJd 3 kG 2 (k) 
J d 3 kG 2 (k) (£(k) + E(-k)) 



G 



Jd 3 kG 2 (k) 

Substituting the explicit form of the spectrum ( |38| ) into (E(k)) G and supposing that 
| g(k) | 2 ~ 1 we obtain: 



/'< 



,(5 

2M 



k 2 



(2nY 



V* XJd 3 kG 2 (k) 



(91) 



Hereupon one can see that if the quantity M is interpreted as an effective mass of the 
excitation then the first term in (|90|) represents the sum of the kinetic energies each of 
the excitations. The last two terms describe the bond energy of the excitations. 

The solution of the equation ([n]) we obtain as an asymptotic series over the coupling 
constant A. Using again the approximate expression for 5 2 , viz. keeping | g(k) | 2 ^ 1 as 



in (|83D we can write equation (|7T|) in the following form: 

XM ' i -(^ |2 



1 



(27T) 



t\m 



k 2 + 5 2 



(92) 



This equation defines 6 2 as a function of A. We will seek solution of ( |92"D in the form of 
the asymptotic series 



5 2 

k 2 + S 2 



ci , c 2 



A ' A 2 



a + aX + — + -^ - 
k 2 a 



X \a + 



X 



c\ c 2 

+ — + — + 
A 2 A 3 



(93) 



By means of the series conversion formulae we can derive 



A 



1 



k 2 + 5 2 a + be 



where e = \, b = a$ + k 2 



ce z 



l--e + 

a 



h l 



a , 



e 2 + 



(94) 



: = ci, ect. Substituting the series expansion Q9"3| ) into the 
equation (]92|), collecting the terms of the same power over A and equating them to zero, we 



will obtain the chain of relations defining the coefficients ao, a, ci, ect. Straightforward 
calculations give the following result: 

M 



a = -(k), o=- v#J 

M /^|fl(k)| 2 (k 2 -(k 2 )) 2 _ _V- 
Cl ~ V* /d 3 A;|(7(k)| 2 " M a l } 



k 2 



;^k 2 u(k) 



2 



/d3fc|^(k)|2 

from which we receive the asymptotic series for 5 2 and //: 

,2 /, 2\ M ^ V * ° 

According to the condition (^t]) 5 > 0. This condition the restricts region of the admissible 
values of M and henceforth Vq. The magnitude of M is determined from the vacuum 
energy Wq minimization condition. 

From the relations ( pop and ( p6[ ) it is easy to see that the nonperturbative and singular 
contributions over the coupling constant A into the energy are defined by the dispersion 
a over momentum distribution of | g(k) | 2 inside the excitation. 

5. Minimization of the Vacuum Energy Wo With Respect to the Transfor- 
mation Parameters. 

We will go on now calculating the vacuum energy Wo and minimizing it with respect to 
the parameters. This procedure corresponds physically to the fact that at the spontaneous 
transition, described by the transformations fllQ|) , the system itself fixes the magnitudes 
of the parameters, choosing the states with the lowest energy. 

As it follows from the relations (p3|), the vacuum energy Wo is equal to: 

Wo = M J d^xe-^eMe^ + 2v^ (M + _A_) + cons t. (97) 

where const, is defined by result of bringing the conterterm AH (2) to the normal ordering 
form with respect to the Heisenberg fields, vq and a(x) are the transformation parameters. 
When we have been obtaining the spectrum of one-particle excitations, as a necessary 
condition of their existence we have found that a(x) = kox . We retain a(x) in the 
expression for the vacuum energy Wq as an arbitrary function determining it later from 
the conditions of the absolute minimum existing for Wo over the variables vq and a(x). 
As it will be demonstrated in a moment we will have the same result, what points to the 
strict correlation between the existence of minimum and the one-particle excitation. 

To calculate const, in ( |9?| ) it is convenient to pass in the relation ( |5^ ) from the fields 
4> a (x) to ip a (x) in accordance with the transformations fllPl) . As a result we obtain: 

-AH{2) = u (ul-vl)e Q ^(x)Us7) + ^)v(x)^(x) + 
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A 



+ u (ul - v 2 ) v(x)e a/3 ij} a (x) e(y) + tztti ^A x ) ~ 



2uy ^(x) e( V ) + 



+ 2u 2 v(x)ij a (x) e(y) 



2V* 



A 



2V* 



if) a (x) + 



v ( x )i'a( X ) 



(9£ 



We can see that the only term in 
the last addendum. Therefore 



const. 



) which does not have the normal ordering form is 
A 



2u 2 / d 6 xv(x) <tl; a (x), e(y) + 



V* 



d xe 



-ia(x) 



2V* 

V* 



e(v)e ia(x) - 4— M 5^ 




(99) 



Together with the representation (|9~T| ) this leads to the following expression for the energy 

-W = p(vo)F(a) + 7(v )A, 



V 



where 



F{a) = i f d 3 xe- ia(x) e(\7 
p(vo) 



Ja(x) 



A 



(k 2 ) + A 



2m ' 2V* 1 
2v 2 (l - 2(1 - t; 2 ) 2 ) , j(v ) = -2t; 2 (l - 2^ 2 ) 



(100) 



The functional F(a) turns out to be essentially positive quantity. Varying it with respect 
to a(x) and putting the first variation to be equal to zero, it is easy to find that an 
extremum of F(a) (minimum in fact) is achieved on the functions satisfying the equation 
\/ 2 a{x) = 0. Evidently, the extremum would define a minimal value of the energy Wq 
provided that p(fo) > 0. Therefore, the minimization of Wo over Vq should be achieved 
within the region where this condition is fulfilled. From the equation for a(x) one can 
easily write out the only solution as a(x) = k x. Let us study now the behavior of Wo 
as a function of parameter v 2 = —k. There are three possible cases corresponding to the 

> O. 



three different regions of -^ '■ 1 + -ttH < 0: 1 + -tt-t = 0: 1 + -,, 

to e(ko) e(k ) ' e(k ) ' e(k ) 

The value of Wq at the minimum point should be negative and, besides, in this point 
p(k.o) > 0. It is seen from the explicit form of p(k) (relation (|100|) ) that this condition is 

-0.3. 



realized within to regions, when k 



> and when —0.5 < k < — 1 + 



The left boundary of this inequality is determined by the condition (^). The equation 
on k is obtained from the condition that the derivative of Wo with respect to k is equal 
to zero and has the form 




0. 



(101) 



Hence we have 



K, 






(102) 
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At 1 + jy^r — > oo,k — > — | ~ —0.375; and at n = — 1 + 4= from ( |102| ) follows that 
1 + ^v ~ 2. Thus, the vacuum energy VFo has the absolute minimum with respect to 
the variables a(x) and Vq, provided that 1 + -^^ > 2, and position of this minimum lies 
within the region 0.3 < Vq < 0.375. The excitation mass M = x ^i is larger then the 



fermionic mass m. There are no solutions inside the region < 1 + jt^- < 2. Finally, 
if 1 + ^£-r < then minimum W is situated at the region of negative values of v%, and 
transformation (plj) corresponds to the hyperbolic rotation. Using expression for A and 



([T00|) we get 

' < (103) 




1m 2V* 
Let us rewrite this relation on the form either 

k% A /(k 2 ) A 



or 



2m + F + tr-5F <0 (104) 



< (105) 




2m 



Hence it follows that 2£'(0) < 0. Thus, the state A + \ 0) will be the lowest state with 
the energy equal to fi. The excitation mass M in this case will be less than m. Further, 
contribution to the bound state energy \x corresponding to the kinetic energy of the pair 
moving as a whole (second term in the sum fl9l|) ) will have positive sign but the third term, 
describing interaction of the excitations - negative sign, in agreement with the physical 
picture. 



Conclusion. 

Finally we want to emphasize some aspects of this paper. The basic and most essential 
one is postulating the dynamical mapping, but the transformations (§) and fllPf ) have been 
used just as a motivation sample to find the mapping (|20|) . 

The next important step is connected with the introduction of the distribution function 
g(k), which allows to define the product of the physical fields at the same space-time 
point. This fact corresponds to the "cut-off" at large momenta, and the free parameter 
V* (spatial volume of excitation), appearing in our model, is related with the "cut-off" 
parameter A cut by the expression: 

4vr VA cut 

Such a smooth way of regularization leads to the fact that the anticommutation relations 
for (j) a (x) (|14D cease to be local, that results, in its turn, in the arising of abnormal nonzero 
anticommutators between the Heisenberg fields. 

As a nontrivial cosequence of " g" -regularization we have the fact that in all dynamical 
characteristics - energy spectrum, vacuum energy density, bound state mass,- the coupling 
constant A is renormalized, such that the effective coupling constant turns out to be y?. 
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The method of dynamical mappings in the form used here can be easily generalized 
on the relativistic model (NJL), at least in separable approximation. This problem will 
be considered in the forthcoming paper. 
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Appendix. 

Here we present the formulas for some physical quantities. 
I. Condensates. 

ul 



(0 | ^ a {x)^{x) | 0) = ^S ap] 

(o | 4>+(x)Mx) I o) = ^:<W; (A.1) 
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II. Casual Green functions for the Heisenberg fields. 



| TM*)4>t(v) I °) = u ° G «k x -V) + ^oV ko(x - y) GS(* - y); 
| T1>+(x)Mv) I 0) = v 2 e-^~^G^(x - y) + u 2 G^(x - y); 



(A 2) 
| Ttf(x)i/>$(y) | 0) = uovoe^, (e" ik ° x G« (x - y) - e-^G%(x - y)) ; 

(0 | T,J> a {x)My) I °) = -^v e a p, (f°*G%{x -y)- e^Gf (x - y)) , 
where G^{x — y), i — 1, 2 are the casual Green functions for the physical fields: 

G%{x-y) = (0 | T<t> a {x)^{y) | o) ; 

G^(x-y) = (0\T<f>+{x)My)\0). (A ' 3) 



III. Superpropagators. 

Fourier image of the function G^l(x-y), by the conventional terminology (see e.g., [[II 



defines the superpropagator presented by the infinite seria of free field propagators with 
initial spectrum e(k). Let us find this representation. We have, for example: 

G ^ (X yj " (2tt) 4 7 cu - £(k) + ze 6 ^ (A.4) 

Extracting the initial spectrum e(k) from the spectrum E(k): 

E(k) = |^(k)| 2 ^ e (k)- Vo 2 e (k-k )-A^^ 

~ u 2 e(k) - v 2 e(k - ko) - A^ = e(k) - E(k), ( A ' 5 ) 



where E(k) = v 2 fe(k) + e(k — k ) + y^j, we come to the expression for the super- 
propagator: 

-i oo 

^(-i) m - 1 (s(k)) m " 1 A m (k), 

1 (A.6) 



w - E(k) + ze m=1 



A(k) 



u — e(k) + ze 



We would like to mention that self-energy part of S(k) is not an analytical function 
over the four fermionic coupling constant, this fact is related with non-analyticity of Vq(X). 
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